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ABSTRACT 


The evolution of weak disturbances in inert binary mixtures is deter- 
mined for the one-dimensional pis* rr problem. The interaction of the dis- 
si‘ Ttive and nonlinear mechanisms is described by Burgers' equation. The 
binary-mixture diffusion mechanisms enter as an additive term in an effec- 
tive diffusivity. Results for the impulsive motion of a piston moving into 
an ambient medium and the sinusoidally oscillating piston are used to illus- 
trate the results and elucidate the incorrect behavior pertaining to the 
associated linear theory. 
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1. INTRODUCTION 


Combined dissipative effects of viscosity, thermal conduction, and 
mass diffusion play fundamental roles in the propagation of disturbancas 
within mixtures of gases. In many practical applications the behavior 
of weak disturbances is of primary concern. Whereas the acoustic, or 
linear, limit leads' to great mathematical simplications, there are situ- 
ations in which nonlinear cumulative effects cannot be ignored. In this 
paper we wish to descrjl‘~ the interaction of nonlinear and the above- 
mentioned dissipative mechanisms as they pertain to the evol ;t’.on of 
weak plane waves in inert binary mixtures. This is done by means of the 
classical piston problem. 

Some of the basics of acoustics for binary mixtures, especially 
with regards to sound absorption, were set forth by Kohler (1949). The 
propagation of spherical acoustic disturbances in binary mixtures, with 
thermal diffusion ignored, was studied by Rasmussen and Frair (1976) . 

For binary mixtures, the study of nonlinear behavior appears to be lim- 
ited to steady-state shock waves. Dyakov (1954) studied weak shocks and 
Sherman (1960) studied both weak and strong shocks. On the other hamd, 
for pure gases, when mass diffusion does not play a role, combined non- 
linear and dissipative effects on the evolution of finite disturbances 
have been the subject of investigation for some time. The classic early 
work is that of Lighthill (1956) , and a summary of other related work 
dealing witn Burgers' equation is given by Benton and Platzmann (1972) . 

A more recent work is that of Halabicky and Sirovich (1973). Also, for 
pure gases, Shidlovsky (1975, 1977) has investigated the evolution of 
shock waves and other regions of tioiiun iformity by means of boundary- 
layer singular-perturbation methods. 


General model equations have been proposed by Blythe (1969) and 
Ockendon and Spence (1969) for the evolution of waves in relaxing gases. 
The thrust of their work was toward inviscid flows, but their models 
could also deal with the viscous-type relaxation under consideration in 
the present work. 

In this investigation, we describe how weak, one-dimensional distur- 
bances in binary inert mixtures evolve with time. In particular, we pro- 
ceed from the basic linear theory and then illustrate how the smallest 
nonlinear terms modify the evolutionary description. For tlie piston 
problem, tne wave front is shown to be governed by Burgers' equation. 

The contributions of the binary-mixture diffusion mechanisms enter as an 

additive term in an effective diffusivity that is a combination of the 
Schmidt number and pressure and thermal diffusion coefficients. Besides 

the mass average velocity perturbation, results are obtained for the spe- 
cies mass-fraction perturbation, which is pertinent to binary mixtures. 
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2. BASIC EQUATIONS 


The equations of motion for mass, species, momentum, and energy for 
an inert binary mixture are (with body forces neglected) 


^ + p div V = 0 


( 1 ) 


DC _ 
^ Dt " 


- div 1 


( 2 ) 


DV ■*-*■ 

p — = - grad p + div T 
Dt 


( 3 ) 


- Dh Dp , j . 

P — = — + T ; e - div q 
Dt Dt 


(4) 


where p, p, h, and V are the density, pressure, specific enthalpy, ar.d 
velocity of the mixture. The mass fractions of the two inert species of 
the binary mixture are denoted by c^ = c and c^ = 1 - c- The rate of 

• 4 — ► 

strain tensor, e, and viscous stress tensor, T, are given by 


If '* t 

£ = — grad V + (grad V) 


(5) 


T = 2ue + X(div V) I 


( 6 ) 


where p and X are the first and second coefficients of viscosity. The 
heat-flux vector, q, and binary diffusion-flux vector, i, are determined 
by means of kinetic theory (Hirschf elder et al., 1954) or principles ot 
continuum mechanics (Landau and Lifshitz, 1959); 

q = - k grad T + jh^ - h^ + ^ 


(In T) 


( 8 ) 


Here k is the thermal conductivity, binary diffusion coeffi- 
cient, and given by 



( 9 ) 


where T is the temperature, and are the specific gas constants for 
species 1 and 2, and and k^ are the pressure-diffusion and thermal- 
diffusion coefficients, given by 



( 10 ) 


k^ E ac(l - c) 


( 11 ) 


The thermal diffusion factor a is usually positive when = c refers to 
the heavier molecular species, but may be slightly negative for excep- 
tional gas pairs. For Maxwellian interaction potentials, a vanishes. 

The set of equations becomes complete with the addition of thermal 
and caloric equations of state. For a mixture of thermally perfect gases, 
we have 

p = Pt|(r^ - R2)c + Rj] , h = (hj - h^jc + h^ (12) 
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3. LINEARIZED EQUATIONS 


Conbider a uniform ambient state denoted by the subscript naught, 
We consider perturbations about this state and write 

P “ Po(l + P') » P - Pp(l + p') , T = T^(l + t')] 


c = c + c 
o 


V/a. 


(13) 


where denotes the frozen speed of sound in the ambient state, 


‘f =YPo/Po 
o 


Y - 1 = 
Y 


(^1 ~ ^2)^0 ^2 

/C - C \c + c 


(14) 


and where C 


and C 


are 


heats of species 1 and°2. 


T, 2 uid position vector, r. 


the ainbient-state constant-pressure specific 
Further, we introduce a nondimensional time, 
as 


I 


2 

t/v 
f o 
o 


-► 

r 


'■‘f ^/^o 
o 


ab) 


such that V H 3/9r (the barred space variables are dimensional). Char.ic- 
teristic Prandtl and Schmidt numbers are defined ns 


P 

r 


p V c /k 

OOP 

^o 


o 


S 

c 


/D, 

o' 12 


(16) 


where 


U 

o 






o ' 


- c 


10 


P. 


20 ' 


+ C 

I* 


20 


(17) 


and K is the hulk nHjdiilua of viscosity in the ambient medium, 
o 

W I ow consider flows that are irrotational (such as occur with 
planar, cylindrical, or spherical symmetry) and introduce a velocity 
potential such that v* = Viji. The linearized momentum equation can thus 
be integrated once, and the governing linearized system of equations 
becomes 


+ V^(j> * 0 , D c' = k y^p' + k„ v^t' 

91 s p^ 


(18) 


p' = Y[y^<{) - If] , p' - p' + t' + N^c' 


M k._ P 


, - # o T r r 

D t' ^ p ^ ^ 2_ v^c' 

P *Y r 9 t S 

c L 


+ k V^p' + k V^T 

P T 

o o 


1 


(19) 


( 20 ) 


where 


M 


o 


S M. 


12 


/I- 


T 

o 



( 21 ) 


N 

o 



Rj - R^ 




D 

s 


c <)t 



( 22 ) 


(23) 


The above equations can be manipulated so that the following single 
seventh-order equation for ({i is obtained: 

Y|1 + I vS.| t v'’|v‘’((. - 


4 V' 


3i«,r - ^2^^ ^ 



TT 


0 


(24) 


where 


a = P 


(1 - I . Y(1 Y sj 


YN k 

op 


+ YP IK 


k_ 


M k, 

o T 


e, E a + S P - Y/1 + N k„ 
1 c r \ on 

= S + P /l + M k^ 1 
2 c r\ oT^I 


( 25 ) 


When k = 0, 
T 


the above linear equations reduce to those of P-asmussen 


o 

and Frair (1976) who studied spherical waves by means of these equations. 

The solution of the linearized equations for a one-dimensional [;is- 
ton moving into an ambient medium requires three boundary conditions at 
the piston surface; (1) the mass mean velocity, (J)^, is equal to th.- pL: - 
ton velocity, (2) the diffusion-flux ncrmal to the i>iston vanishes, 

1 • n = 0, and (^) cither the temperature nr tlie nomal lieat flux, q • ;*■ , 
is Rfiecifiofl. Th(? solution is obt<iin»-<i analogously to the r.]>heric.il ex- 
plosion problem of Rasmussen and I'rair (l')7b). Tin; solution can )>• 
written as a part that describes .i diffused wave front and two coupletl 
parts that describe heat and mass-diffusion boundary layers adjacent to 
tilt piston face. The );eat and ma.ss-di f fusion boundary layers spread from 
the piston face as the square root of time whereas the position of the 
bulk v^ave front travels like the lime itself. Wiien we fix our .ittention 
on the wave front for large times, therefore, v^e c.ui ignore the boundary 
layers since they are greatly oiitdistancerl by the wave Iront. Tie* a:;ymp- 
totic long-time wave behavior can be obtained by I.apl ace-t r inaf rirn teoh- 
niqiies (Itasmussen 1975). For th*' velocity and mass-f rac;t ion ju-r t ur b it ion 
Wo have far the impuln.ive motion of a pi .ton, wher.- 'I'^fx - 0, t > 0) 


U /a, , 

|) f o 


^ v' (x,T) 

erfc 

b(x - T) 

4b^x _ 

+ e erfc 

b(x + T) 


2 a. 
fo 


/T 





c'(x,T) 


“phpo " ''I - » vl 


fo c 


bx 

»4r^ 


exp 


2 2 
b'‘(x - T)'' 


( 26 ) 


( 27 ) 


where 


b = 


(^) 

\ r c ' 


( 28 ) 


The above results show that the wave front is centered at x = T and 

spreads out with time. The peak of the mass-fraction perturbation at 

X = T dies out like x”** and thus ultimately vanishes. If we define the 

thickness of the wave front as 6 = | — —I /l9v^/3x| , we obtain for larq 

\a / * 'max' 


times 


a- 6 

fo 


/nr 

b 


( 29 ) 


Thus the wave front in the linear theory spreads out like the square- 
root of time. It is known, however, that a constant piston speed should 
generate a steady-state shock front. The continual spreading of the wave 
front and the dying out of the mass-fraction perturbation are incorrect 
behaviors that arise out of the omission of the convective nonlinearities. 
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4 . NONLINEAR INTERACTION 


In order to delineate the balance between the dissipative transport 
terms and the nonlinear terms, let us first renormalize the variables ap- 
pearing in (24). We introduce a characteristic length L and define non- 
dimensional time and distance as 

T = a^ t/L = ET , r = r/L = Er , (30) 

o 

where the parameter E = L is to be regarded as small. We note that 

o 

V = eV- In the new variables, nth-order terms are proportional to E^, 
and we can rewrite (24) and display the lowest two orders as 

‘ - MttL 


Thus as the viscous (transport) effects go to zero, the classical wave 
operator prevails, and the first-order correction for small E arises be- 
cause of the next higher-order derivatives. The terms of order E in 
(31) correspond to the seventh- and six-order derivatives in (24) . 

Equation (31) does not account for the nonlinearities in the prob- 
lem. We account for tlie lower-order nonlinearities by replacing the per- 
turbation quantities in (13) by p* = Ep, and so on for the other primed 
variables, and utilizing the normalizations (30). The full nonlinear 
equations then lead to the following equation for the perturbation po- 
tential ; 


TT 


= e 


TT 


.vi. 

P s 
r c 




* |(Y - * o(e^) 

TT 


( 32 ) 


In Equation (32) the second-degree nonlinear terms appear to order e, as 


I 

2 

shown, and the third and higher degree nonlinearities appear to order e . 
The nonlinear terms of order e in (32) are the same as occur in inviscid 
potential theory. Thus, the first-order corrections to inviscid acoustics 
are the sum of the first-order linear viscous correction and the first- 
order nonlinear inviscid correction, which of course might have been an- 
ticipated a priori. Our goal now is to examine Equation (32) with the 
2 

terms of order e ignored. 



5. PISTON PROBLEM 


The boundary condition for the one-dimensional impulsive motion of 

a piston is (x = 0,T) = (U /a^ )f(T), where U is the characteristic 

p fo p 

speed of the piston and f(x) is of order unity. The characteristic 

length is taken to be L = V /U , and it follows from (30) that 

o p 

T = U a- t/v„, X = U x/V , and e = U /a-: . The small parameter £ is 
pf o p o p f 

o ^ ^ o 

thus the frozen Mach number of the piston, and we note that v* = 

A series solution of Equation (32) by a straight-forward expansion 
in powers of E will lead to secular behavior such that the first-order 
correction becomes as large as the zeroth-order term when x = o(c ^). 

We account for this behavior by introducing new variables defined as 

C=x-T , X = ex . (33) 

We further expand ({) in a series of the form 

4)(x,T;e) = 4 >q(C»X) + E4>j^(CfX) + . . . (34) 


Expressions (33) and (34) applied to Equation (32) lead to the following 
equation for <})^(5,X), which has been integrated twice with respect to ^ 
cuid the functions of integration set equal to zero: 



+ — (p (p 



2P S 


r c 



(35) 


We now note that Ur* is the lowest order velocity contribution, 
o o^ 

that is. 


V 


/ 


eu 

o 



(36) 
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If we make the further definitions 


n 


= T 1 „ ^ Y + 1 
2 ^ 2 


eic 


V* 



(y + 1)P s 

r c 




(37) 


then Equation (35) can be ejq)ressed in the standard form for Burgers’ 
equation ; 


U + U U 
o o o 

n 


v*u 


(38) 


The nonlinear terms of order e in (32) lead to the convective nonlinear 

term U U in (38), and the dissipative terms of order e in (32) lead to 
o o^ 

the diffusion term U in (38) . The origin of the diffusivity factor 
V* can be seen by identifying the factors and original 

linear Equation (24) . To lowest order, the other perturbation variables 
are found to be 


p' = EYUq + o(e^) 

p' = eu^ + o(e^) 

t' = e(Y - i)u^ + 0 (c^) 

■ - S K * - ““t I"o, * 

c I o o ] C 


(39) 
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5.1 Impulsive Compression 


For impulsive motion of the piston into the fluid, the function 


U must satisfy the boundary condition U (C,0) = where H(Q 

o o 

is the Heaviside unit step function. The solution to (38) in tliis case 
is derived in Chapter 4 of Whitham (1974) and is 

ertc 

' /4n* ' 


Uc,(5*»n*) - 




(40) 


where = ^/V* and r|* = H/V*. The binary-mixture contribution appears 
in the parameter V* and hence amounts to a stretching of the coordinates. 
The derivative U is given by 


u (^*,n*) 


_ Jr I!1 


2 n*) 

(C*-n*) /4n* i„^^, 

2 o 




(41) 


The function U^((^*,r|*) = “ T*,c*x*) is plotted in Figure 1 as 

a function of x* = x/v* for various times T* = T/V* and for 
c* = (y + l)e/2 = 0.1, where n* = C*x*. In the early stages, evolution 
of the wave front is dominated by diffusion and hence flattens until a 
balance is achieved by the steepening effects of nonlinear convection. 

The ^-derivative of the velocity, U , which is a measure of the slope 
of the wave front for small e, is snown in Figure 2. This figure indi- 
cates how the maximum slope decreases and the wave front broadens as 
time increases. The evolution is dramatic for T* £ 14 shown in Figures 1 
and 2, but Tor larger times changes occur more slowly. The maximum value 


of -U approaches 0.125 as T* > ", but at T* - 140 the maximum value of 
Or* 



-U has decreased to only 0.141, for e* = 0.1. Ultimately, diffusion 
and nonlinear convection balance and a steady-state shape is attained 
with -U having a symmetric shape about a maximum value of 1/8. 

The function U also describes the mass-fraction perturbation, c , 
as seen from expressions (39). For the heavier species |Yk + (Y “ 1 

I o| 

is negative, and hence c* is negative in the wave front- Thus, the mass 
fraction of the heavy species decreases in the wave front, and the mass 
fraction of the light species increases. The maximum value of these spe- 
cies perturbations is skewed toward the piston side of the wave front. As 
time increases, the distribution broadens and ultimately becomes symmetric, 
achieving a steady-state form. 


5.2 Thickness of Steady-State Wave Front 

As T* “, the wave front evolves into a steady-state form. The 
thickness 6 of the wave front can be defined as 


6 = 


1 



max 


( 42 ) 


For T* and E -*■ 0, we obtain 


V ° (^1 - ^2 

(Y + UP S 

V r c 


8 

Y - 1 + P YN 

r o / 


Y + 1 

P Ski 

P Y 

-T 


r c p ' 

o 


o ' 


( 43 ) 


Note here that P and S are defined in terms of the reduced viscosity 
r c 

V = (2u + X ) /p . Formula (43) is identical to previous results 

o o o o 
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obtained for weak steady shocks, such as by Sherman (1960) when rewritten 
in our notation, which is a partial check on the correctness of our pres- 
ent analysis. Note that the parauneter V* is related to the parameter b 

2 

in the linear theory. Equation (28), by the relation 2(y + l)b V* = 1. 
Thus linear theory yields the correct combination of terms arising from 
viscous, thermal, and mass-diffusion dissipation that contribute to the 
breadth of the wave front, but raised to the wrong power. 

The combination N k can be written as 

° Po 


N k 

o 




(44) 


where m^^ and m^ are the molecular masses of the species. Since this com- 
bination is always positive, it can be seen from (43) that pressure and 
thermal diffusion always tend to broaden the wave front, the more so the 
greater the disparity in masses of the species. As pointed out by Slicrman 

(1960), k and k^ are usually of opposite signs and tend to counteract 
^o o 

one another. 


5.3 Oscillating Piston 


For the oscillating piston problem, the boundary condition is 

v*(0,T) = e sin u)T, where £ = U /a^ as before, and ui is the nondimens ional 

frequency defined in terms of the physical frequency, (J, by E v (o/a^ U . 

o f p 
o 

The initial condition on the Burgers' function is thus ^^(£[,0 = 0) = 

- sin to^. By means of Fourier series expansions, we find the solution of 


(38) to be (Benton and Platzmann, 1972) 


1 ) 

o 


(C,n) 


4(OV*^^nA^exp|-n^co''V*nlsin (ncofj 


n=l 


I 2'2 I ~ 

1 + 2^Aj^oxp{-n V*n)f-Os (n(o£,) 


n=I 


(4'.) 
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where 


A (Cv*) H 
n n 

and t is the modified Bessel function of the first kind of order n. The 
n 

appearance of the factor V* delineates the bxnary -mixture, or mass diffu- 
sion, effects. 

~2 

When u) V*n » 1, the term n = 1 dominates the series, and thus in 
the limit of ultimate decay expression (45) becomes 

~ ~ -2 - 

U (CfH) 4u)V*A ( 0 JV*)exp(-O) v*n)sin (47) 

o 1 

Aside from the exponential damping, the amplitude of the disturbance in 
the far field is a function of the frequency of oscillation. This is a 
nonlinear effect euid is akin somewhat to the dependency of the amplitude 
on frequency in the nonlinear oscillations of spring-mass systems. The 
amplitude function 4(I5v*A^ (alv*) is plotted as a function of wv* in Figure 3. 
The function is a maximum of unity when WV* -> and decreases as tou* be- 
comes smaller. When uiv* » 1, expression (47) approaches the form 

~2 

U^(5,n) ''' - exp(-u) U)^ (48) 

This is the solution that is obtained from the linear theory, described 
by Equation (27), in the limit eiv* « 1. Thus, the linear theory ac- 
cording to Equation (27) is valid in the final stage of decay when 
0)U* » 1 and e << (WV*) Except for thir. special limit, nonlinear ef- 

fects always make a contribution, even in the final stages of decay. 

Even though small, the nonlinear effects are cumulative. At high 


( 2wv* ^ ^ oilin '^* ) 


(46) 
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frequencies the steepening of the compression part of the wave tends to 
c^mcel the flattening of the expansion part, whereas for low frequencies 
the compressions and expansions are less closely spaced emd hence lead to 
a nonvcuilshing cumulative effect in the far field. 

The evolution of the wave disturbances just before final decay Ccin 
be exeunined by keeping the term n = 1 in the denominator as well as the 


numerator of (45) and ignoring all terms for n ^ 2. Then we have, for 

to^v*n ~ 1, 

-2 

exp (-U) V*ri) sin 

— T : 

1 - — 1 exp (-W V*n) cos 



The wave forms produced by expression (49) are shown in Figure 4. The 

~2 

curves are shown for (i) v*n = 1 and wv* = “, 0.5, 0.25. The linear sine 
wave is represented by tou* = The values uiv* = 0.5 and 0.25 represent 
the effects of the nonlinearities, which yield wave forms that are 
flatter in the middle and steeper on the ends. As n for oiv* fixed, 

all the curves evolve toward the sine-wave shape. 

The other perturbation variables are also given by expressions (39) , 
the pressure, dt;.jity, and temperature being proportional to the velocity. 
The mass-fraction perturbation, c* , is proportional to the ^'derivative 
of For the behavior approaching ultimate decay, v*n ~ 1, the de- 

rivative of (49) yields a valid description: 



o 
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The behavior of expression (50), and thus of c', is shown in Figure 5 for 

u V*n = 1 and 0)V* = 0.5, and 0.25. The limiting value loV* * « yields 

a damped cosine wave that corresponds to the associated linear theory. 

The perturbation c' is thus 90 degrees out of phase with U . When uv* is 

o 

finite, the magnitude of U , and thus c', is larger in the comoression 

part of the wave than it is in the expansion part. This is a nonlinear 

effect. If c' is identified with the heavier species, such that 

(Yk + (y - 1) k,j. } is negative, then we find that c' is negative in the 
*^o o 

compression part of the wave and positive with a smaller magnitude in the 
expansion part of the wave. It also can be seen that more time is spent 
in the expansion part of the wave than in the compression part. As n 
for U)V* fixed, all the curves tend toward the cosine-wave shape. 
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6. CONCLUDfNG RKMAFiKS 


Tlio intoraction of nonlinear atul dissipative nechanisms for w-*ak di.j 
turbanoos in inert binary mixtures lias been inve.il iqati.-d . A.s for pure 
gases, the basic description of one-cumensional problems is by means of 
Burgers' equation. The effects of the binary mixture an- del iru?atifd b*/ 
the etfective ditfusion coefficient V* which contains the effects of 
Schmidt number and pressure and thermal diffusion coefficients associated 
with binary mixtures as well as the I’randtl number and ratio of specific 
heats associated also with pure gases. The binary-mixture mechanisms 
tend to broaden the wave fronts and to modify the amplitude-frequency re- 
lation and damping rate of oscillatory disturbances. The rnas; -f ract ion 
perturbation is skewed toward the jiisfon side of evolvin*) slu, k front:;, 
but becomes ;;/mmetr ica I as the sti'ady- .tate i;; approached . For osi.iila- 
t.ory motions, tin- magnitude of th" ract ion p'Tt ;irbal i on is laroer 

in the compression part of the wave and of malli-r mac-.itude and op[«site 
sign in the expansion part of the wave. 

The nonzero mass-fraction perturbations suggest Unit prt-ssure and 
ther. i diffusion may have significant effects in chemically reacting 
mixtures where nonequil \br i urn effects play essential roles. Research 
towards these ends may yield significant results for the behavior of det- 
on.ition waves and combustion noise. 

Perturbation methods applied to hi' ^r-or ’er iorm:i of t m uliuear 
Equation { i2) m.iy yield mote gi'neral moi. . e luatii li hin .duraets’ <' ;ua- 
tion. Also spherical and cyiindric.il d ist urF i.cct n 'ht bi* studii'd 

The;..- 


fruitfully by tfii;; ai.proach. 


aie ar.»as tor future ri';;e,irch. 



uirc 1. Kvolution of 




Figure 3. Variation of Amplitude 
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Fiqure S. Mass-Fraction Perturbation Wave Forms for Final Decay 
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